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Abstract 

The pp — ► {pp} s 7T differential cross section has been measured with the ANKE 
spectrometer at COSY-Jiilich for seven proton beam energies T p between 0.5 and 
1.97 GeV. By selecting proton pairs with an excitation energy of less than 3 MeV it 
is ensured that the final {pp} s system is in the state. In the measured region of 
®pp ~ 18°, the data reveal a forward dip for T p < 1.4 GeV whereas a forward peaking 
is seen at 1.97 GeV. The energy dependence of the forward cross section shows a 
broad peak in the 0.6-0.8 GeV region, probably associated with A(1232) excitation, 
and a minimum at 1.4 GeV. Some of these features are similar to those observed 
for the spin-isospin partner reaction, pp — > dir + . However, the ratio of the forward 
differential cross sections of the two reactions shows a significant suppression of 
single pion production associated with a spin-singlet final nucleon pair. 
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Single pion production in nucleon-nucleon collisions, NN — > NNir, is one of 
the principal tools used in the investigation of NN dynamics at intermediate 
energies [1-3]. Because of the large momentum transfers involved, even close to 
threshold, such a meson production process is sensitive to the short-distance 
part of the NN interaction. 

The pp — > dir + reaction has been the subject of extensive experimental study 
with the measurement of many spin observables, as well as of the differential 
cross section from threshold up to several GeV. However, the information that 
this provides is restricted to final NN states with spin S = 1 and isospin T = 0. 
On the other hand, the pp — * {pp} s 7r° process is kinematically very similar to 
this provided that the excitation energy in the final proton pair is very small. 
In this case, due the the Pauli principle, the protons must be in the singlet 
1 <S'o state, i.e. have quantum numbers (J p , T, L) = (0 + , 1,0) compared to the 
(1 + , 0,0 and 2) of the deuteron. Pion production in the two cases therefore 
involves different transition matrix elements so that a combined study of the 
two processes should yield greater insight into the reaction dynamics. 

If we consider the pp — > dn + reaction as the limit of triplet {pn} t production in 
the pp — > {pn} t 7T + channel, where the strong final state interaction produces 
the deuteron [4], the ratio of 7r° to tt + cross sections will provide information 
on the relative strength of spin-singlet to spin-triplet production, i.e. give 
information on the relative probability of pion production at short distances 
in channels with different spin orientation of the final nucleons. Because of the 
smallness of the signal, attempts to identify spin-singlet production directly 
from data on the pp — > pnn + cross section have only yielded upper limits [5-7]. 

A small value of the singlet-triplet ratio is expected near threshold since s- 
wave isovector pion rescattering is absent for pp — > ppn° and heavy (cu) meson 
exchange provides the largest driving term [8]. A small value of the ratio is 
also predicted for energies around 0.4-0.6 GeV since the >S-wave A(1232)iV 
intermediate state that dominates the pp — > dn + cross section is forbidden by 
conservation laws in the pp — > {pp} s 7T° case. The theoretical situation at higher 
energies is largely open. The position is rather similar on the experimental side 
since, away from the low energy domain, the only other published data in the 
^So conditions were limited to energies T p < 0.425 GeV [9]. 

We have previously reported a measurement of the pp — > {pp} s 7i° differential 
cross section obtained using the ANKE spectrometer [10] for the single beam 
energy of 0.8 GeV [11] . Here the two protons were detected at small angles 
with respect to the incident beam and cuts were made such that the excitation 
energy E pp in the final pp system was less than 3 MeV. Under these conditions 
we expect the final {pp} s pair to be almost purely in the singlet ^So state. It 
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was found for this beam energy that the angular variation was rather strong 
but that at all the measured angles the cross section was orders of magni- 
tude smaller than that for pp — > dir + . To study the energy dependence of 
these effects, we present here further measurements taken over a wide range 
of energies, 0.508, 0.625, 0.700, 1.100, 1.400, and 1.970 GeV. 

The magnetic spectrometer ANKE is placed at an internal beam station of 
the COSY cooler synchrotron of the Forschungszentrum Julich. Fast charged 
particles, resulting from the interaction of the proton beam with the hydro- 
gen cluster-jet target [12] and passing through the analysing magnetic field, 
were recorded in the forward detector (FD) system [13]. The FD system in- 
cludes multiwire proportional chambers for tracking and a scintillation counter 
hodoscope for energy loss and timing measurements. The triggers employed re- 
quired the crossing of the two planes of the hodoscope by at least one charged 
particle (SP-trigger) or by two particles (DP-trigger) [14]. When the DP- 
trigger was used for data taking, the SP-trigger ran (prescaled) in parallel 
for luminosity measurement and calibration purposes. The tracking system 
provided momentum resolution o p jp pa 1% in the range of interest and res- 
olution in the excitation energy of a(E pp ) pa 0.2 — 0.8 MeV for events with 
E pp < 3 MeV. 

Additional details of the experimental setup and the measurement procedure 
are to be found in refs. [11, 13, 15, 16]. The data at 0.8 GeV and above were 
taken during a single beam-time run, whereas the lower energy results were 
obtained from other ANKE calibration runs. 

From measurements of the momenta of two charged particles in ANKE, the 
pp — ► ppir° channel was isolated by determining the missing mass M x in the 
reaction. In more than 80% of cases where two fast particles were detected, 
the tracks passed through different counters of the forward hodoscope. For 
these events the particles could be clearly identified as protons on the basis of 
the timing information. The difference of the arrival times of the two particles 
measured by the counters was compared with the time-of-flight difference de- 
duced from the measured momenta assuming that the particles both had the 
mass of the proton. For the remaining pa 20% of events, the hypothesis was 
made that the two particles were indeed protons. As already shown for the 
0.8 GeV data [11], the missing-mass distributions for both classes of events are 
very similar with only a slightly enhanced background when no timing infor- 
mation was available. The two sets were therefore combined in the subsequent 
analysis. 

Figured] presents examples of the measured missing-mass distributions in the 
7T° region for events with E pp < 3 MeV. In addition to the 7r° peak, a rise of 
counts is seen on the right hand side due to two-pion production which seems 
to be largest at 1.1 and 1.4 GeV. This background gives little contribution in 
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Ml [(GeV/c 2 ) 2 ] 

Fig. 1. Missing-mass-squared distributions of the pp — > ppX reaction for M x in 
the 7T° region. The curves show the fits to the experimental spectra in terms of a 
Gaussian plus a straight line. The data near M% = were excluded from the fit 
since, as indicated here by the 0.625 GeV results, there may be some pp — ► {pp\ s l 
events in this region. 

the 7T° region and the data were fitted as a sum of a Gaussian and straight line. 
The region close to = was excluded from the fit since, as discussed for 
the 0.8 GeV data [11], there is the possibility here of some contribution from 
the pp — > {pp}s1 reaction. In all cases the peak position was consistent with 
m^o to within the experimental uncertainties of about ±10MeV/c 2 . Events 
within ±2er of the central value were retained for the determination of the 
pp — > pp7r° differential cross section. The numbers of 7r° events deduced in this 
way are given for the different energies in Table [TJ 

The luminosities recorded in Table Q] were obtained by measuring in parallel 
elastic proton-proton scattering using the SP trigger. The ANKE setup detects 
the fast proton produced by this reaction for cms angles between about 10° 
and 30°. However, to avoid regions where the acceptance changes rapidly with 
angle, only data from the range 15° < < 24° were retained for 1.4 GeV 
and below while at 1.97 GeV the interval 17° < 9™ < 27° was selected. After 
corrections for acceptance, the numbers of events in several bins of 8p m were 
compared with the values of the elastic differential cross section taken from the 
SP07 solution provided by the SAID phase shift analysis [17]. The precision of 
these predictions was checked by looking at experimental data at small angles 
from which it was assessed to be typically about ±4%, though a little larger 
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Table 1 

Summary of experimental conditions: L- m t is the integrated luminosity at the beam 
energy T p and N^ g /N^o is the fractional background under the peak, where N n o is 
the number of pp — > {pp} s ir° events registered in the angular region 6™ < 15°. 



at 2 GeV. 

A full Monte Carlo simulation of the ANKE spectrometer has been developed 
within the framework of the GEANT program [18]. This allowed us to estimate 
the acceptance factors for different bins of the kinematic variables and hence, 
on the basis of the luminosities given in Table [lj to evaluate differential cross 
sections. 

The experimental distributions in the final proton-proton variables for our 
data at all energies are very similar to those published at 0.8 GeV [11]. The 
excitation energy E pp spectrum is well described in terms of an S"-wave pp final 
state interaction [4], provided that the Coulomb force is included. For E pp < 
3 MeV the angular distribution of the pp relative momentum vector in the 
final pp rest frame is consistent with isotropy, as expected for the production 
of a ^So pair. 

Although the coverage in the diproton angle 6™ with respect to the beam 
direction is rather limited in the ANKE spectrometer, we already noted that 
the data at 0.8 GeV showed a strong dependence on this angle [11]. The same is 
true for the other energies shown in Fig. [2] where, since the two initial protons 
are identical, the data are plotted as functions of cos 2 9 pp ". For energies of 
1.4 GeV and below, the results show a forward dip whereas at the highest 
energy the cross section is maximal in the forward direction. This is perhaps 
an indication that the reaction mechanism changes with energy as different 
intermediate nucleon isobars are excited. 
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Fig. 2. Differential cross section for the pp — > {pp} s 7r° reaction for proton-proton 
excitation energies E pp < 3MeV. The results at different beam energies are shown 
in terms of cos 2 The straight lines represent fits according to Eq. (pQ) with the 
resulting parameters being given for all energies in Table [2j 

The angular distributions have been fitted with the linear form 



= (T + ffi sin 2 9™ 

d <>;•'/ 



and the values obtained for the forward differential cross section o"o and the 
slope parameter o\ are recorded in Table [2j 

Figure [3] summarises the results from this and a previous experiment [9] by 
showing the energy dependence of the forward differential cross section and 
the slope parameter o\. Although the statistical error on the 0.508 GeV point 
is very large, the pp — > {pp} s ^° cross section data suggest a significant rise 
from this energy to a maximum in the 0.6-0.7 GeV region, where A(1232) 
production is expected to be particularly strong, followed by a monotonic 
decline. However, the 2 GeV point lies much above this trend and this, together 
with the slope information shown in the lower panel, indicates that the data 
have entered here a different domain. A similar tendency is seen in the pp — > 
dn + forward cross sections, also shown on the figure, where the cross section 
displays a broad minimum around 1.4 GeV after leaving the region influenced 
by the A(1232) isobar. However, 7r + production rises at lower energies due 
to contributions from the excitation of the A(1232)iV system in a relative S 
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Fig. 3. Upper panel: Energy dependence of the forward differential cross section 
for the pp — » {pp} s tt reaction with E pp < 3MeV. The closed circles represent 
the results from the present experiment while the triangles show the low energy 
CELSIUS data [9]. For comparison we show also the corresponding cross section 
for the pp — > dir + reaction. For energies up to 1.4 GeV this is represented by the 
dashed line taken from the SAID parameterisation [20] whereas at higher energies 
the Akemoto et al. data are shown as open circles [21]. Lower panel: Slope parameter 
o\ of the pp — > {pp} s ir data, as defined by Eq. (pQ). For both panels the solid curve 
corresponds to the predictions from the A^iV/A(1232)A r model of Niskanen [19]. 



wave, which is excluded on spin-parity grounds in the 7r° case. 



The behaviour of the slope parameter a\ shown in the lower panel of Fig. [3] is 
rather different, with a much smoother variation and a change of sign between 
1.40 and 1.97 GeV. It is interesting to note that data on the pp — > dn + differ- 
ential cross section also show a small forward dip in the 0.6-1.4 GeV energy 
range [20] and, furthermore, that the sign of the forward slope changes in the 
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Table 2 

The forward differential cross section o~o and slope parameter o\ of Eq. ([I]) for the 
pp — > {pp} s 7T reaction with i?pp < 3 MeV. The results at 0.425 GeV and below are 
from Ref. [9]. Also shown are the forward cross sections (177,-+) for pp — > dir + [20,21] 
and the ratio i?(-7r°/7r + ) of the forward pp — > {pp} s vr° topp — > {pn}t7r + cross sections 
integrated up to 3 MeV excitation energy. 

1.9-2.1 GeV region [21,22]. 

The only theoretical estimates of the cross section and slope parameter for 
the pp — > {pp}s^° reaction for the small E pp kinematical domain have been 
made by Niskanen [19]. The cross section predictions shown in Fig. [3], which re- 
flect the delicate interferences between contributions from NN and A(1232)A r 
intermediate states, do illustrate a displacement of the A(1232)iV peak up- 
wards by about 100 MeV compared to pp — > d,7i + . Although this agrees with 
our observations, Niskanen's predicted minimum around 0.7 GeV is in com- 
plete contradiction to our findings of a maximum in this region. Furthermore 
we see no sign of any big peak in the slope <j\ around 0.6 GeV so that any 
agreement with our 0.8 GeV data seems fortuitous. 

To facilitate the comparison of tt + and 7r° production, the values of both 
the pp — > {pp} s tt° and pp — > dir + forward cross sections are given in Ta- 
ble |2j However, to compare the strengths of pion production leading to the 
S-wave spin-triplet and singlet AW states, we require rather data on the 
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PP ~ * {pn}t7r + channel integrated up to 3 MeV excitation energy. Using final- 
state-interaction theory, this can be approximated for low pn excitation ener- 
gies by the pp — > dir + data through [23] 



^(pp ^ {pn} t TT + )xi^(pp ^ dn + ) x J 



da 



k 2 



Tra t (k 2 + 0$ 



dk 



« 0.096 x ^-(w -> rf7T+) , (2) 
dll 

where k is the relative momentum in the final NN system, k^^Jm^ = 3 MeV, 
and 0% = ttinB, with B being the deuteron binding energy. Though there are 
some deviations from this, due in part to the pn tensor force, the approxima- 
tion gives a plausible description of experimental data [7] . 




T p [GeV] 

Fig. 4. Energy dependence of the ratio i?(-7r° /vr + ) of the forward pp — > {pp} s 7r° and 
PP ~^ {fwi}t7r + cross sections integrated up to 3 MeV excitation energy. The curve 
is drawn to guide the eye. 

A common feature of the R(ir°/'K + ) ratio seen in Fig. H] is its relatively low 
value. Below 1.4 GeV it is typically a few per cent, which is significantly smaller 
than the trivial spin-statistics factor of 1/3. These conclusions are not changed 
significantly if the difference between the spin-singlet and triplet final state 
interactions are taken into account [24]. 

The dip in the ratio in the 0.5 GeV region is in part a reflection of the fact 
that the S"-wave A(1232)iV intermediate state is forbidden for ir° production 
whereas it plays a vital role in the case of the 7r + . However it is not at all clear 
why R(tc /tc + ) is small in the 1 GeV region where the contribution from the 
A(1232) maximum is much reduced. On the other hand, the much larger value 
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at 1.97 GeV is yet another indication that the reaction mechanisms might be 
different at higher energies. 

In summary, we have measured the differential cross section for the pp — > 
{pp} s Tr° reaction at seven beam energies from 0.5 to 1.97 GeV under the spe- 
cific kinematic conditions where the proton-proton excitation energy is below 
3 MeV and the cms angle between the diproton momentum and the beam axis 
is less than 18°. The observed form of the E pp spectra and isotropy of the an- 
gular distribution in the diproton cms are consistent with the assumption that 
the two final protons are in the ^So state. Except for the highest energy, the 
data all show a sharp minimum in the forward direction and the ratio there of 
the differential cross sections for the pp — > {pp} s ^° and pp — > dn + reactions is 
below 1%. The situation changes radically at 2 GeV when a forward maximum 
is observed with a much enhanced value of i?(7r /7r + ). This must be a reflec- 
tion of a different reaction dynamics at high energies. A broad peak around 
0.6-0.7 GeV is observed in the energy dependence of the zero-degree differ- 
ential cross section. The data are consistent with the predicted displacement 
of the A(1232) maximum to higher energies for 7r° production but the other 
features of our results are not explained by a microscopic model [19]. This 
could, for example, be due to phase differences between intermediate NN and 
A(1232)iV contributions and much theoretical work still remains to be done. 

We have attempted to study the relative amplitudes for the production of 
spin-singlet and spin-triplet final NN states by using a simplistic model that 
links the pp — > dn + and pp —> {pn} t ft + cross sections. The resulting i?(7r°/7r + ) 
ratio shows that singlet production in the forward direction remains small at 
energies even above the A(1232) excitation region. 

Further data are clearly required in the 0.5-0.6 GeV region where the Niskanen 
model should be more reliable than at higher energies to see if there is indeed 
any trace of his predicted structure there. The transition region above 1.4 GeV 
is also of interest. Does the spin-singlet suppression continue when higher 
isobars enter? 

For a final state there are only two spin-amplitudes for the pp — > {pp} s ^° 
reaction. Since we observe a very strong angular variation over our acceptance 
region, it is then likely that the analysing power will also be significant there. 
It is therefore planned that in the near future measurements of A y will be 
carried out at COSY [25] and eventually spin-correlation studies might be 
achievable [26]. This ensemble of data should provide a significant challenge 
for theory in this, one of the simplest pion-production reactions. 
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